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ABSTRACT 

We present a study of X-ray ionization of magnetohydrodynamic (MHD) 
accretion-disk winds in an effort to constrain the physics underlying the highly- 
ionized ultra-fast outflows (UFOs) inferred by X-ray absorbers often detected in 
various sub-classes of Seyfert active galactic nuclei (AGNs). Our primary focus is 
to show that magnetically-driven outflows are indeed physically plausible candi¬ 
dates for the observed outflows accounting for the AGN absorption properties of 
the present X-ray spectroscopic observations. Employing a stratihed MHD wind 
launched across the entire AGN accretion disk, we calculate its X-ray ionization 
and the ensuing X-ray absorption line spectra. Assuming an appropriate ionizing 
AGN spectrum, we apply our MHD winds to model the absorption features in 
an AMM-Aewton/EPIG spectrum of the narrow-line Seyfert, PG 1211-1-143. We 
hnd, through identifying the detected features with Fe Ka transitions, that the 
absorber has a characteristic ionization parameter of log(^c[erg cm s“^]) ~ 5 — 6 
and a column density on the order of Nh — 10^^ cm“^, outflowing at a charac¬ 
teristic velocity of Vc/c ~ 0.1 — 0.2 (where c is the speed of light). The best-£t 
model favors its radial location at ~ 200Ro {Ro is the black hole innermost 
stable circular orbit), with an inner wind truncation radius at Rt — 30i?o- The 
overall K-shell feature in the data is suggested to be dominated by Fe xxv with 
very little contribution from Fe XXVI and weakly-ionized iron, which is in a good 
agreement with a series of earlier analysis of the UFOs in various AGNs including 
PG 1211+143. 


Subject headings: accretion, accretion disks — galaxies; Seyfert — methods: numerical 
— galaxies: individual (PG 1211+143) — X-rays; galaxies 
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1. Introduction 


Blueshifted absorption lines are among the most common spectral features seen in 
the spectra of accreting compact objects across a large dynamic range in black hole (BH) 
mass from the upermassive BHs of active galactic nuclei (AGNs) to stellar-mass black 
holes of galactic binary systems. In the former case, approximately 50% of Seyferts and 
quasars (QSOs) exhibit absorption signatures in the UV band (e.g. ICrenshaw et al.l 1199911 
with a similar fracti on 50%) of Seyfer t I’s showing Huesh ifted absorption features in 
their X-ray spectra (IRevnolds et al.lll997l: iGeorge et al.lll998h indicative of an underlying 
physical link between these two outflow components. A small fraction (~ 10%) of the 
radio-quiet QSOs further shows substantially blueshifted UV resonance lines, referred to 
as broad-absorption-lines (BALs); these are mainly Giv/ Mgll (high/low-ionizat ion) at 
velocities of n/c ~ 0.04 — 0.1 where c is speed of light (e.g. Grenshaw et al.l2003[l . 


X-ray spectroscopy plays a fundamental role in the study of AGN absorber properties, 
because, in distinction to the UV transitions, the X-ray ones span a much wider range in 
ionization parameter ^ (i.e. the ratio of photon to electron fluxes); thus, within the span of 
1.5 decades in photon energy (~ 0.3 — 10 keV) one can sample atomic transitions that cover 
5 decades in ^ (e.g. from neutral to Fexxvi) and presumably a large range in physical 
length scales. Typically, the so-called X-ray warm absorbers (WAs) have characteristic 
local columns of Nh ^ 10^^ cm~^ and ionization parameter in the range — 1 ^ log^ ^ 4 at 
LoS velocities of v/c < 0.01 (e.g. [Reynolds fc Fabian! 1199511 presumably originating from 
sub-pc to pc scales. A rich spectral diversity in soft X-ray regime (< 2 — 3 keV) with a 
large number of ion ic transitions affor ds the statistical studies of their X-ray absorption 
line properties (e.g. iBehar et al.ll2003h . Among them, the absorption measure distribution 
(AMD) can be used as a global measure of the density of the radiation-absorbing gas 
along a line-of-sight (LoS). AMD is the differential hydrogen-equivalent column Nh per 
decade of i.e. (iNn/dlog^, and it is computed from the measur ed columns of a variety o f 


ionic species of several elements spanning a large range in £ (e .g. ISteenbrugge et al.l 12005 


Holczer. Behar fc Kaspill2007l: lBeharll2009l:lDetmers et al.ll201111 . The AMD determination 
in a number of radio-quiet Seyferts seems to indicate, to zeroth order, a similar global 
column distribution (i.e. roughly a constant AMD), implying a wind density n(r) that is 


similar in all of them a nd decreases like n{r) oc r ^ with radius r (e.g. iDetmers et al.ll2011 


Holczer &: Beharll2012[l . 


Furthermore, in recent years, exhaustive X-ray studies of Fe K-shell transitions in 
AGNs by XMM-Newton and Suzaku have revealed the presence of another outflowing 
component in the Seyfert spectra, typically identified as highly-ionized high-Z ions such as 
Fexxv/Fexxvi with H-equivalent columns of Nh ^ 10^^ cm“^, high ionization parameter 
(log'C ^ 4) at near-relativistic outflow speeds of n/c > 0.03, named for this reason ultra-fast 
outflows (UFOs). The detected UFOs appear to be ubiquitous across both radio-quiet 
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Seyferts like PG 1211+143 (see alsolReeyes et al.ll2009l: iPounds et al.ll2003l: iPounds fc Page 


20061: iTombesi et al.l l2010ai 1201 lal. l2012al: iGofford et al.l I2013L l2014i) and radio-loud ones 


(e.g. 3C 111 and 3C 120, 3C 390.3, 3C 445) wi th a likely association of their properties to 
the radio spectra flTombesi et al.ll2010bl l2011bl) . The higher X-ray content and increased 
wind ionization of the latter suggests that strong X-ray photoionization apparently does 
not inhibit the launch of such fast winds. A detailed study of the properties of X-ray 
absorbers in a samp le of 23 AGN using high resolution X-ray spectroscopy was conducted by 
Blustin et al.l fl2005[) . to conclude that most of the X-ray absorbing matter is launched from 
large radii (the AGN molecular torus) with kinetic luminosities that are only a small fraction 
of the AGN budget. In addition to these Seyferts and nearby QSOs, optically/UV-bright 
broad-absorption-line (BAL) QSOs and their variants (e.g. non-BAL and mini-BAL QSOs) 
apparently show similar X-ray UFOs but wit h even higher velocities, up to vjc ^ 0 .7 — 0.8 
in extreme case^ such as APM 08279-1-5255 f Ghartas et al. 2002, 2003, 2007 . 2009 ). 


PG 1211-1-143 is a bright quasar at redshift ^ = 0.0809 flMarziani et al.lll996l ) with 
X-ray luminosity ~ 10^^ erg s“^ in 2 — 10 keV band for Hq = 75 km s“^ Mpc“^ and Galactic 
hydrogen-equivalent column density Nu = 2.85 x 10^° cm“^ (Murphy et al. 1996). It is 
an optically bright quasar with a prominent “Big Blue Bump” that results in a relatively 
steep optical/UV-to-X-ray flux ratic@ (oox = —1.45). Among other spectral features, the 
first AMM-Aetwfon/EPIG/RGS observation of PG 1211-1-143 in 2001 revealed the first 
evidence of a highly-ionized UFO with mass flux and kinetic energ y comparable to tha t 


of the accretion mass r ate and bolometric luminosity, respectively fjPonnds et ahl 12003 


Pounds fc Reevesll2009l) . although we note that ot hers have reached different conclusion s 
depending on how to tre at the baseline continua flKaspi fc Beharl 120061: iGallo fc Fabian 
20131 : IZoghbi et al.N2015l ). Their analyses detected several strong absorption features 
identified as blueshifted Ka transitions of G, N, O, Ne, Mg, S and Fe. The properties 
of the Fe feature in particular, imply a wind column density of Nu ~ 5 x 10^^ cm“^ 
at velocity vjc ^ 0.08 and ionization parameter log^ ~ 3.4. A second observation of 


PG 1 211-1-143 with AMM-Ae wfon/EPIG/RGS in 2004 (IPounds fc ReevesI 120071 ) and in 
2007 fjPounds fc ReevesI l2009l) have again detected similar UFOs, implying their persistent 


^While the winds of the typical high-velocity UV transition (La, Giv etc.) BALs may 
be driven by line radiation-pressure, it is a challenge for this scenario to accelerate the 
highly-ionized, near-relativistic X-ray UFOs. 

^The spectral index pqx = 0.384 log(/ 2 kev// 25 oo) measures the X-ray-to-UV relative 

o 

brightness where / 2 kev and /2500 are, respectively, 2 keV and 2500 A flux densities (Tanan- 
baum et al. 1979). 
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presence despite its highly variable X-ray spectra. A more detailed spectr al analysis of su ch 


outflows has been recently performed using the XMM-Newton data (e.g. iTombesi et al 


2010al. l2011al: lPoundsll2014l) to confirm their presence, in agreement with the earlier results. 
Final ly, a more recent observation with Suzaku/XlS has also re vealed the same UFOs as 
well ( Reeves et al. 2008 : Patrick et al. 2012 : Gofford et al. 2013 1. 


Despite the long-known UV/X-ray WAs combined with an increasing number of 
statistically-signiflcant detections of the X-ray UFOs, the detailed geometrical structure of 
these ionized winds, including formation and acceleration processes, are poorly constrained 
to date. Yet, each of these issues is crucial to the comprehensive picture of accretion-powered 
phenomena in accretion/outflow physi cs. Plausible launching mechanisms for general 
outflows include radiation-driven (e.g. Proga. Stone, fc Kallmani 2000 : Proea fc Kallman 


Fukumura et al. 


Fnknmura et al. 


2005 


2004: iNomura et al.ll2013L i n the context of UV BALs in lu minous QSOsl, thermally- driven 


(e.g. 

Begelman e 

al. 

1983 

) and magnetically-driven 

(e.g. 

Blandford & Pavne 

1982: 

Konigl & Kartiel 

1994 

. Gontoponlos & Lovelace 

1994 

, hereafter GL94: 

Ferreira 

1997; 


2010al. hereafter FKCBlOa: iFukumura et al.ll2010bl. hereafter F KCBlOb; 
2014. hereafter F141. There have a l so been hyb rid models (e.g. 


Everett 


Ohsuga et al.l 120091 : lOhsnga fc Mineshigell201ll : IProga 1 20031 ) that attempt to explain 


an AGN phenomenology associated with inflow and outflovu. With increasingly improved 
fully-numerical schemes, various extensive simulations have been made in the context o f 


me aisK-wma 

Forth & Fendt 

jcenari 

2916: 

o lor ID magn 

Mnrnhv et al. 

erican 

201 Q: 

V ariven le.g 

Stnte et al. 1 

. rnc 

mii 

niz er ai. 

Stepanovs 

/uun: renar /uun: 

& Fendtll2ni4h and 

(hi radiation driven (e.g. 

Proga & Kallman |2004: 

Nomura et al. 

2013: 

Higginbottom et al. 


2014: iHagino et al.ll2015l) . Although the acceleration mechanism(s) of the observed winds 
remains uncertain, the magnetic origin seems to be favored over the radiation pressure 
one according to the latest time-dependent hydrodynamic simulatio ns coupled with 


multi -dimension Monte Garlo calcula tions for radiative transfer (e.g.. iHigginbottom et al. 


2005 


Everett 


2005 


2014 . but also see IHagino et al.ll2015r) and from UV/X-ray observations (e.g. iKraemer et al. 


Grenshaw fc Kraemerl l2007h 


Miller et al. 2006 


This may also be the case for Galactic 
King et al.ll2012l. 120141) . One should note that certain 


binaries (e.g. 

phenomenological outflow models, with emphasis on individual spectral features, such 
as the Fe K-shell transitions, are able to reproduce the properties of certain prominent 
transitions such as their EW and their LoS velocity (e.g. ISimll2005l : ISim et al.ll2008l . 12010 : 
Tatum et al.ll2012h without providing, however, a global dynamic wind perspective. 


To the best of our knowledge, none of the existing wind models, whether semi-analytic 
or numerical, have so far been able to deliver a practical prescription for the observed X-ray 


^The derived values of large Nh and v of certain X-ray UFOs presumably originating 
from smaller radii are a serious challenge against line-driven and thermally-driven scenarios. 
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absorption features; i.e. local properties (i.e. column, ionization state, velocity) of the WAs 
and UFOs together with a global picture of the outflow physics (i.e. density/ionization 
structure from smaller scale to larger scales and geometrical properties as a whole). From 
methodological viewpoint, most models £t the properties of specihc features, i.e. column 
and velocity, implementing xspec/xstar to obtain the ionization parameter and the 
velocity of the plasma associated with specific transitions, with little concern on how these 
£t within a global model of the AGN outflows. 


The spirit of our recent works (i.e. CL94; FKCBlOa; FKCBlOb; iKazanas et ahl 12012 : 
F14) has been exactly the opposite, in that we begin with a global MHD wind model and 
use the X-ray spectroscopic observations to determine the global properties of these winds. 
In this paper we employ a similar philosophy in an attempt to model the observed Fexxv 
UFO in PG 1211-1-143 within the context of the well-dehned MHD-driven wind models 
referred to above. Our study allows us to constrain explicitly some of the defining MHD 
wind parameters in the spirit of a model-driven approach. Our deeper goal is a better 
understanding of the underlying physical structure of the observed winds from a global 
standpoint. Within this framework, WA and UFO features are generically identihed as 
belonging to the same wind structure that spans the entire domain of the AGN accretion 
disk. We briefly describe the essence of the MHD-driven winds in §2 along with our 
methodology for constructing a grid of simulated line spectra for subsequent data analysis. 
In §3 we show our preliminary results based on a 60-ks XMM-Newton/'EPIC spectrum of 
PG 1211-1-143 deriving the best-£t values for the primary model variables. We summarize 
and discuss the implications of the model in §4. 


2. Ultra-fast Outflows in Stratified MHD Disk-Winds 

2.1. The Magnetized Disk-Wind Structnre 

Following FKGBlOa and FKGBlOb for the computational prescription of magnetically- 
driven disk-wind models under steady-state, axisymmetric conditions, we seek new insight 
into their structure from the observational data. We apply our model assuming the observed 
X-ray UFO signatures in AGNs (i.e. Fexxv/Fexxvi resonance transitions)0 are produced 
by X-ray photoionization of MHD winds launched off of an accretion disk. The detailed 
characteristics of the model discussed in FKGBlOa and FKGBlOb will be briefly described 


“^The model, however, is not restricted to Fe K-shell transitions and can be extended in 
general to include other ionic features detected in AGNs and black hole binaries. 
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here. Geometric and physical properties of the wind in the model are primarily governed 
by two conserved qnantities along a wind streamline, namely the particle-to-magnetic flnx 
ratio Fo and angnlar momentnm Ho- The former one, Fq, determines predominantly the 
wind kinematics and the latter, i7o, generally dictates the global wind structnre in the 
poloidal plane. The fnndamental qnantity of axisymmetric MHD is the magnetic stream 
fnnction T(r, 6*), assnmed to have a self-similar form '^{r,6) = with Tq the 

poloidal magnetic flnx through the fiducial innermost disk radius at R = Rq. ^(0) is its 
angular dependence to be solved for and g is a free parameter that determines the radial 
dependence of the poloidal current. The scalings of the poloidal magnetic stream function 
carry over to the rest of the wind properties of which we show only the magnetic field, 
velocity and density 


B{r,e) = {R/Roy-^B{e)Bo , (1) 

v{r,9) = {R/Ro)~^^‘^v{9)vo , (2) 

n{r,9) = {R/Rof'^~^h{9)B‘^^v-‘^m-\ (3) 

with the momentum-balance equation 

p(v-V)v = -Vp - pV$g -F ^(J X B) , (4) 


where rup is the proton mass and p is plasma mass density. 

The dimensionless angular functions denoted by tilde must be obtained from the 
conservation equations and the solution of the Grad-Shafranov equation (the force balance 
equation in the 6*—direction) with initial values on the disk (denoted by the subscript “o”) 
at {R = Ro, 9 = 90°). The density normalization at {Ro, 90°), setting h(90°) = 1, is given in 
terms of the dimensionless mass-accretion rate rha (normalized to the Eddington accretion 
rate Me = LeI(?-, see FKGBlOa) by 


n 


O 


O'tRs 


(5) 


where <7^ is the Thomson cross-section, is the ratio of the outflow rate in the wind to rha 
and Ro is assumed to be on the order of the Schwarzschild radius Rs- The Thomson depth 
T{rha) of the plasma at the innermost disk radius is further scaled by the dimensionless 
mass-accretion rate rha with normalization as T{rha) = rhaTo, which leads to 


Un = 


fwFla,o 

O'tRs 


= 5 


fwFia,c 


X 10^^ cm ^ 


( 6 ) 


where we have introduced an effective mass-accretion rate rha^o = idiaTo as it is difficult to 
decouple one from the other from observations. In this p aper we consider one of the hducial 
wind solutions, model (A), from iFnknmnra et af.l (120141) . as a baseline wind model, by 








choosing q = 0.93 (i.e. n oc fw = ^ and Tq = 10 representing an optically-thick disk 

of rha o = 10 at i? = Rs (low case r denotes the radial distance in 3-space, while R along 
the disk surface) . Here, we o nly highlight the esse nce of the model; details can be found 
elsewhere (CL94; FKCBlOa.b: iKazanas et ahll2012l: F14). Formally the self-similar winds 
extend from r = 0 to r —)• cxo, however, physical considerations restrict these to a hnite but 
broad range in r. So we choose the dimensionless factors ft-, /r to denote the inner and 
outer truncation radii of our winds on the disk surface by 


Rt = ftRs , Rt = ItRs ( 7 ) 

where the value of ft is to be constrained by the X-ray data while /t 3> 1, typically ~ 10®. 
Once launched, the asymptotic wind speed in this solution is found to be Vp/vo ~ 4 at 
r/Rs < 10® (see F14 for details). 




Fig. 1.— (a) Input ionizing SED for photoionization calculations for the MOD of fcTbbb = 
10,30,50 and 70 eV with ctox = —1-5 and F = 2. (b) An example of calculated Fexxv 
(hydrogen-equivalent) local column distribution Mh for different sets of {9, /cTbbb) as labeled. 


2.2. Photoionization of Disk-Winds 

With the dimensionless, mass-invariant wind structure for a given rha o and a 
viewing angle 6, the only signihcant difference in the wind ionization properties across 
objects of different luminosity comes from the spectral energy distribution (SED) of the 
accretion-powered luminosity L = rhaLE^ where Le = 1.25 x 10^® Mg erg s“® is the 
Eddington luminosity with Mg the black hole mass in units of 10®Mq and e ~ 0.1 is the 
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accretion effici ency. While in FKCB lOa we nsed a simple power-law spectrnm of the form 
Fj, oc (e.g. Sim et ah! 2008 . 2010 ). here we consider a mnlti-component SED consisting of 


a mnlticolor-disk (MCD) with an innermost temperatnre of fcTbbb and an X-ray power-law 
of photon index E (with a low-energy cnt-o ff at 50 eV and a high-e nergy tnrnover at 200 
keV) normalized to the MCD by aox (e.g. lEverettI l2005l: ISimI 120051) . a more appropriate 
SED for bright Seyferts snch as PG 1211-1-143. The ionizing Inminosity (X-ray plus EUV) 
is then Lion — O.IL ~ 1.25M8 x 10^^ e rg s~^ for a relative ly high accretion rate of nta = 1 as 


suggested in the earlier analyses (e.g. iPounds et ahl 120031) 


For a characteristic Seyfert SED, we set T = 2 (see 

Fig. 5 in 

Tombesi et ah 

2011 a 

for a homogeneous sample of 42 radio-quiet AGNs and 

Pounds & Reeves 

2009) and 


aox = —1.5 (adopted from NED and 

Blustin et ah 20051) while leaving an inclination angle 

9 and the disk temperature kTuhh as 

free parameters to be determined bv PG 1211-1-143 

UFO observations fe.g. Pounds et ah 

2003; 

Tombesi et ah 

2011a). It should be noted 

that, in agreement with FKGBlOb, 1 

Blustin et ah f 

2005 

) find that only more negative 


values of aox allow higher velocity absorbers (i.e. Uout ^ 10, 000 km s“^) based on their 
analysis of phenomenological and physical properties of the detected warm absorbers using 
high-resolution X-ray spectroscopy of a sample of Seyfert 1 type AGNs. 

Given the wind density normalization Uq through photoioniza tion balance is 


computed radially outward employing xstar flKallman fc Bautista! 1200 ll. v2.2.1bnl3) by 


setting the SED of Fig. [T^ as the ionizing spectrum at the innermost radius; the radiation 
transport in the wind is done by discretizing the the radial wind coordinate using a large 
number of cells in radial direction for a given angle 9 (typically with Ar/r ~ 0.1 allowing 
to treat each radial cell as a plane, yielding 50-70 radial zones; see FKGBlOa). We apply 
xstar in the first zone to compute the ionization equilibrium of the plasma and also its 
opacity and emissivity. Then we use the output of this zone as input for the next one and 


the absorption sp ectra with the Voigt function 
Hanke et ahl 120091 ) defined as 


T.g. 


Mihalas 

1978; 

Kotani et ah 

2000 


H{a, u) = — 


o-y 


dy 


TT 7-00 (w - yf + 0? 


( 8 ) 


where we use a = r£;/(47rAi/£)) with being the Einstein coefficient and Ai//) the line 
Doppler broadening factor. The dimensionless frequency spread about the transition 
frequency Vo is given by m = (z/ — i'q)/Ai'£). Note that, in order to compute the flux in 
lines whose thermal width is narrower than the computational frequency grid (especially in 
cases of multiple lines within a given frequency spacing), the parameter vturb (typically 
~ 1000 km/s) of xstar is employed to provide line width consistent with the produced 
flux over the grid size. However, our wind model provides, instead, a well-defined velocity 
shear AV with a corresponding radial velocity difference Avd between two adjacent radial 
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cells; we employ this velocity instead of vturb to define an eqnivalent Aud = c)Avd-i 

a value consistent with the underlying wind kinematics (see FKCBlOa for a detailed 
numerical prescription). 

Using the ionic column over a radial cell of width Ar as a function of 

ionization parameter ^(r;6*) obtained with xstar under ionization and heating-cooling 
balance, we can compute the wind opacity T,^{r,d) of any given photon energy and at any 
given point with wind velocity v{r]6) from the relation 

Ty{r, 9) = crphoto,i/(u 0)Nion{r, 9) , (9) 

where the line photoabsorption cross section Uphoto,!/ at frequency z/ is given by 

o-photo,;. = 0.001495 cm^ , (10) 

Aud 

and /ij is the oscillator strength of the transition between the i-th and j-th levels of an ionic 
species. Finally, we construct a two-dimensional grid of baseline spectra for 9 G [30°, 70°] 
and fcTbbb e [lOeV, 70eV] for density normalizatiorfl Uo = 5.1 X 10^^ cm ^ (rhao = 10). 
Here we introduce the quantity Mh dehned as the number density of Fe xxv ions divided 
by the Fe abundance and multiplied by the width of our local radial grid size Ar. Some 
of the calculated Mh (assuming solar abundances) for four sets of 9 and fcTbbb are shown 
as a function of the wind velocity v/c in Figure [T)d. Considering this hgure, it is reminded 
that the velocity decreases with increasing distance r and decreasing ionization parameter 
^ for a given LoS angle 9. The reader should note that Mh does not depend monotonically 
on velocity because, at small r (and high n), a good fraction of Fe is fully ionized, while 
at larger distances (i.e. low velocities) the Fe ionization drops precipitously. Finally, the 
total Mh, NniFexxv), is found by integrating Mh over r along a given LoS. As seen, for 
a given ma,o, the normalization of the LoS column depends primarily on the inclination 9 
while the location of peak Mh (he. where Mh is maximum), for a given aox and Uo, is 
mainly determined by the disk temperature fcTbbb- Such a correlation is also discussed in 
FKCBlOb. 


2.3. A Spectral Grid of mhdwind Models 

For the model described in §2.1-2.2 and with the density normalization hxed at 
77-0 = 5 X 10^^ cm“^ at 9 = 90° (i.e. rha^o = 10) we construct a grid of mhdwind models 

®Note that the wind density n(r, 9) oc f{r)g{9)no has both radial /(r) and angular g{9) 
dependencies. g{9) (Fig. 2a of FKCBlOa) decreases by factor of ~ 10"^ for 0° < 6* < 90°. 
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Fig. 2.— Radial profiles for photoionized wind temperature log(T[K]) (on the left ordinate) 
and velocity vjc (on the right ordinate) as a function of ionization parameter ^ for Q = 50° 
and fcTbbb = 30 eV. 

by varying the three parameters {0, kT\^hh, ft)- By solving the radiative transfer along 
radial rays with the wind photoionization as described above we first obtain the column 
distribution for Fexxv and Fexxvi in the ionized wind as a function of ionization 
parameter ^ (r, 6 ) or equivalently as a function of the wind velocity v along a given LoS 
angle 6 (see FKCB10a,b) as shown in Figure [Tb. A different aspect of the wind is shown 
in Figure [2] where we demonstrate the wind velocity n(.^) and temperature T(0 under 
ionization balance for 6 = 50° and fcTbbb = 30 eV. The wind velocity in this framework 
scales like v oc as discussed in FKCBlOa, FKCBlOb and F14. A sample grid of the 
simulated Fexxv absorption line profiles is shown in Figure Eb for various combinations of 
{9, fcTbbb, Rt) within the range considered here. As seen, the dependence on each parameter 
can be probed by the spectral shape in terms of the depth of trough and line shift. A 
comparison between Fexxv and Fexxvi is shown in Figure Eb where we explore the 
spectral variations for different truncation radius ft = Rt/Ro by setting fcTbbb = 40 eV and 
6 = 50°. One should note that the small fluctuations in the model line prohles in Figure El 
is not real and caused by the size of the radial bin we choose in our model. 
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Fig. 3.— (a) Simulated Fexxv absorption features illustrating the diagnostic capability of 
our model for 6 = 30°, 40°, 50° and 60° (from top to bottom) with thin curves denoting an 
input thermal disk spectrum with kTm = 30 eV and thick curves being 70 eV. The dashed 
curve is given by /cTbbb = 30 eV but for Rt/Ro = 10. (b) A similar K-shell feature due to 
Fe XXV (black) and Fe xxvi (red) for various truncation radii Rt as listed in Table [H 

3. Preliminary Comparison with the PG 1211+143 Data 


3.1. XMM-NewtonfEPlC Data 


We use an XMM-Newton spectrum of PG 1211+143 (obsID: 0112610101 


the E PIC-pn camera for an approximately 60-ks duration on 2001 June 15 flPonnds et al 


obtained with 


20031 ). for which a detailed data reduction procedur e and observed spectral and tempora l 


2006 

: Kaspi & Behar 

2006; 

2009 

: Gallo & Fabian 

2013; 


Tombesi et al. 2011a : Bachev et al. 2009 : Pounds fc Reevesi 


Ponndsll2014j . reference therein). Earlier analyses of the UFOs, 


typically identihed as either Fexxv and/or Fexxvi, seem to imply an estimate on the 
column density of Nh ~ 10^^ cm“^, velocity of n/c ~ 0.1 — 0.15 and ionization parameter 
of log^ rsj 3 — 5, although an alternative view may also be conceivable claiming that the 
observed Fe K absorption feature can be attributed to several consecutive low charge states 


of Fe (see, e.g.. lKaspi fc Beharll2006h . 
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Table 1. Model Grid of mhdwind Component 


Primary Parameter 

Range 

Viewing Angle 6 (degrees) 

BBB Disk Temperature /cTbbb (eV) 

Disk Truncation Radius log/t = log{Rt/Ro) 

30°,40°,50°, 60°, 70° 

10, 30, 50, 70 

0,0.3,0.6, 0.9,1.2,1.5,1.8 


Assume Mg = 1, Oox = —1-5, F = 2, no = 5 x 10^^ cm ^ and Lion = 1-25 x 10^^ erg s 
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3.2. Spectral Modeling for the Fe Ka UFO 


Here we perform a spectral analysis of the UF Os previously detected in the 60 ks 


XXM-Newtonl'EP\C spectr um of PG 1211+143 (e.g. iPounds et ahll2003l: iPounds fc Reeves 


20091 : iTombesi et alJl2011al ). Focusing on the hard X-ray absorption feature identified 


as Fe K-shell resonance transition in the data, we implement our MHD-wind model, 
mhdwind, into xspec as a multiplica t ive tab le model as discussed in §2. We follow the 
analysis procedure in Tombesi et al.l f 2011ail where the 2 — 10 keV band is modeled with 
an underlying continuum power-law (po). To fit the Fe Ka absorber, however, we replace 
the phenomenological xstar component by our MHD-wind model mhdwind. The symbolic 
spectral form reads as “phabs*(po+zga)*mtable{mhdwind}” where we have used the 
previously estimated values of parameters. Galactic absorption due to neutral hydroge n 
column (phabs) Nh = 2.85 x 10^° cm“^ f Murphv et ahl lOOOll . F = 2 ( Tombesi et ah 


2011ah. and black hole mass, M = IO^Mq, based on the earlier estimatets ( Kaspi et ah 


2000l: iBentz et al.ll2009h . Attributing the pronounced emission line at ~ 6.5 keV (in the 


rest-frame) to fiuorescence from the disk, our mhdwind is constrained simultaneously with 
a red-shifted gaussian component zga of XSPEG in which the line width is set to be 
CTpe = 0.15 keV whose exact value has little infiuence on our end results. 


We explore two cases by simultaneously considering both the Fexxv and Fexxvi 
transitions: Model (A), where, the innermost radial extent of the wind at 6 = 90°, is 
equal to Rq — Rs and model (B) where this restriction is relaxed. Figure shows the 
best-fit for each of models (A) and (B) in comparison with the no mhdwind model. We set 
77-0 = 5 X 10^^ cm“^ (n7a,o = 10) while varying 9, fcTbbb and ft. The best-fit values are listed 
in Table 2. We have used different values of Uq in our calculations but its effective role is 
to change slightly the depth of the trough. In the current wind model, both visual and 
statistical inspection favors model (B) by = 2.3 (table 2) in which the wind does not 
originate at the fiducial radius Ro on the disk surface, but at Rt = ftRs with /* > 1. 

In model (A), where ft = 1 is assumed, we obtain our best-fit for values 6 = 40° 
(pegged) and fcTbbb = 30 eV with jv = 200.84/129 with mhdwind which is a statistically 
significant addition to the continuum (with an improvement of = 34.1 for two 
additional parameters). In model (B) we relaxed the restriction on the wind truncation 
radius /* = 1 of model (A). Table [3] shows a list of various characteristic radii in this model. 
Our analysis yields a best-fit model with 6 = 50°, fcTbbb = 38 eV and ft = 10^'"^®, as shown 
in Figure IHi and Table 2, where we obtain = 198.54/128 which is more significant in 
comparison with model (A). We note that the model spectrum now has a sharper edge on 
the bluer side of the feature as required in data. The total column, Mr = 1.2 x 10^^ cm“^ 
from model (B) is comparable to the previous estimate with xstar model although our 
wind is continuous rather than discrete. 


As a measure of assessing the Fexxv absorption wind properties we first calculate a 
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Table 2. Summary of our best-fit mhdwind model parameters for PG 1211-1-143. 


Parameter/Model 

Model (A) 

Model (B) 


Fexxv/Fexxvi 

Fexxv/Fexxvi 

9 [degree] 

40.0 0 

49 8'''3-27 

A:Tbbb [eV] 


QQ -1 +4.55 

^Fe [keV] 

O-o"±-0.080 

'3-^^-0.073 

'^max 

n /nniQ^^'^^23 

u.uyo_g Q]^Q 5 /u.uiy_Q QQ 3 ^ 

U.ZOO_qj 24/ U-Ud^_0.024 

log{rc/Rs)^ 

9 /9 

^•^0-0.161/^•^-*--0.16 

O q7'+0.48 /-I oq+0.54 
^*'^*-0.35/-'-•^^-0.25 

log (^c[erg cm s“^]) 

r (^1+0.149 /r ^(^+0.147 

5.3i;°i3/5.80lHr 

Vc/c ^ 

0.099l°:°23/0.165tH?3 

0.115lO:Oi?/0.208;Hl3 

Mu [cm-2] /1022 # 

4.04l°;f4/5.94l°,s2 

19 1 +5-30 /-I f* ^7+5.72 
-L^‘-L_7.56/ '-8.63 

log (Rt/Rs) 

0 * 

1.48;0;°f 

jv (with mhdwind) 

200.84/129 

198.54/128 

(from phabs* (po+zga)) 

-34.1 

-36.4 


^ The value is pegged. 

^ The characteristic LoS radius Tc where wind Fexxv opacity (see eqn. (9)) is maximum along 
a given LoS angle. 

The characteristic value (“c”) is evaluated at the LoS position r = Vc- 
tt LoS-integrated total Fe xxv column density. 

* The value is fixed. 
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log Nnicm 


23.4 

23.2 

23.0 

22.8 

22.6 


40 


Rest-Frame Energy [keV] 


45 50 55 60 

Inclination Angle 0(deg) 


Fig. 4.— (a) 60 ks XMM-Newton/EPlC-pn spectrum of PG 1211+143 fitted with model 
(A) (blue) and (B) (red) in comparison with no mhdwind component (in green) varying three 
primary parameters (6^, fcTbbb,/t) with Uq = 5 x 10^^ cm“^ (i.e. ma,o = !)• (b) Calculated 
contours for characteristic radius log(rc/i?s') = 1.5, 2, 2.5, 3, 3.5 (solid; upper right to lower 
left) and wind velocity Vc/c = 0.2,0.15,0.1,0.05 (dashed; top to bottom) with the best-fit 
Fexxv colnmn in model (B). Color indicates the net Fe colnmn density A/"^. See Table |2]for 
details. 


characteristic radius at which the wind photoelectric absorption column for the Fe xxv 
transition becomes maximum for a given LoS inclination angle 9. At this radius we compute 
the other physical quantities listed in Table 2. Notice that the total Mh (in nnits of 10^^ 
cm“^) is defined as the local colnmn density integrated over the LoS distance. 

Along the LoS of the valnes of 9 obtained by onr fits (see Table 2), the wind is 
both Thomson thin and also thin at the Fe energies. Because, as argued earlier, the 
Fe XXV /Fe xxvi line opacities are non-monotonic fnnctions of the radial coordinate in these 
directions, we dehne a radius Tc along each of these LoS at which the line(s) opacity (ies) 
is (are) maximized (given by the entry Xmax of Table 2). In fact, these coincide with the 
maxima of Mh of Fig. lb. Because of the smoothness and continnity of Mh with r, ^ or u, 
the absorption of X-ray photons begins at r -C and extends over more than one decade 
in radius. One should hence bear in mind that a given absorption featnre in our models 
does not correspond to specific, nniqne wind component. 

In order to examine in more detail the mnlti-parameter space spanned by {9, kMM) 
we interpolate the wind variables such as velocity Vc and characteristic radius Vf. as shown 
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in Figure Hb where color shows total column log(A/'^[cm“^]) for Fexxv with contours of 
radius logirdRo) (solid), and contours of Vc/c (dashed). The best-£t model (B) for Fexxv 
is indicated by a dark dot. It should be reminded that the best-£t characteristic values 
(i.e. TcVcic) are simply constrained at the most opaque radius (rmax = 'r(r = Tc)) of 
the absorber. The neighboring plasma at r < Tc (i.e. v > Vc) and r > Vc (i.e. v < Vc) 
contributes also progressively to the formation of the absorption feature thus there is no 
single wind velocity nor column density in our model. This is a characteristic feature of the 
continuous wind model which is fundamentally different from a single-component absorber 
model often employed in a phenomenological analysis. The corresponding confidence 
contours for model (B) are shown in Figure O where the primary variables 6, fcTbbb and ft 
are constrained. 




Fig. 5.— Confidence contour plots for (a) temperature fcTbbb and inclination 9 and (b) 
truncation radius Rt and temperature fcTbbb at the confidence level of 68%, 90%, and 99% 
also indicated by the best-£t model (B) (cross). See Table [2] for details. 

In the context of the current model, the best-£t model (B) is spatially identihed 
as illustrated in Figure [6^ where the calculated fiducial wind structure in the vicinity 
of the black hole is superimposed showing the normalized number density n(r, 9) in 
color, the velocity held (white arrows), the magnetic held lines (solid thick), the contours 
for density (dashed) and the Alfven surface (white line). In this simplihed approach a 
geometrically-thin disk is situated in the equatorial plane at 6* = 7r/2. As discussed earlier, 
the faster portion of the modeled Fexxv/Fexxvi absorber (i.e. the bluer side of the trough) 
and the slower one (i.e. the redder one) are respectively located at r < and r > r^ along 
each LoS and they all progressively contribute to produce the observed absorption feature 
(both in depth and width). 
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Fig. 6.— (a) Geometrical identification of the observed Fexxv/Fexxvi from XMM- 
NewtonfEPlC spectrum of PG 1211+143 with mhdwind model (B). We show the global 
magnetic held lines (thick solid), normalized wind density \og{n{r,d)/n^ax) in color and its 
contours (dashed), velocity held (white arrows) along with the inferred location of Fexxv 
and Fe xxvi. See Table [2] for details, (b) A schematic to illustrate the MHD-driven AGN 
winds. 


In terms of the energy budget of the observed UFO, by using the outhow density prohle 
with n oc 1/r (in the present model actually n oc 1/r^’^^) in this wo rk for a supermassive 
black hole mass of M = Mg (IKaspi et ahl l2000l: IPeterson et al.l 120041) , a mass-outhow rate 
associated with the Fe xxv line can be estimated as 


r-200 


Mout (Fe xxv) = 2'Kbm^ 


n{r, 6)vz{r, 9)rdr , 


AirhmpnocR^gVpe xxv 2 : 


1/2 


2.56Mq yr ^ I —- 


exxv-^Fexxv ) 


Tin 


0.4/ 15 X 1011 / VlO^Me 


h 


M 




0.1 




200 


1/2 


(11) 

( 12 ) 

(>3) 


where x = r/Rs and the upper limit of integration is indicative of the distance to Fexxv 
location; i.e. Tc = rppxxv (see tab. 2). This value is consi stent with the earlier estimate 
of ~ SMq yr“i (jPounds et al.ll2003l: IPounds fc Pagell2006l) . Since the corresponding local 


mechanical power is given by 


^(local) ^ J\;^(local)^2 ^^ ^ ^-1/2 


out 


(14) 


the local kinetic power of the Fe xxv outhow is dominated by the inner outhow radius Rt 
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yielding 


Eout(Fexxv) ~ -Mout(FexxvX^^ 


2 X 10^^ erg s 


-1 


(15) 


a val ue comparable to power the observed X-ray luminosity 


1044 gj.g g 1 


flPounds et al. 


2003h also potentially providing a large impact on AGN feedback process at large scales 


(e.g. 


Crenshaw fc Kraemerl 120121) . A similarly la rge outflow power has been made, for 


example, to other bright AGNs such as PDS 456 flReeves et al.l 120031: iNardini et al .112 0151) 


4. Summary &; Discussion 


We have demonstrated, by modeling its XMM-Newton spectrum, that MHD-driven 
winds with n oc r“", a ~ 1, originally proposed to account for the X-ray WAs in 
Seyferts, can also encompass the UFOs, i.e. the high-velocity X-ray absorbers of the bright 
Seyfert PG 1211-1-143. The absorber’s properties of PG 1211-1-143 as manifest by the 
Fexxv/Fexxvi transition properties, are determined mainly by the wind mass flux rha^o, 
the disk temperature /cTbbb and the observer’s viewing angle 9. By producing a grid of 
model K-shell Fe absorption lines, appropriate to photoionized MHD winds, we found 
the that the absorbers’ physical conditions are well constrained by our models. Thus the 
Fe XXV and Fe xxvi properties are respectively given by the location of maximum opacity at 
Vc/Rs = 234 and 66, LoS velocity Vc/c = 0.115 and 0.208, ionization parameter log.^c = 5.31 
and 5.80, total H-equivalent columns Afn = 1.21 x 10^^ cm“^ and 1.67 x 10^^ cm“^, 
with the wind truncated at radius ft = Rt/Rs = 10^'^® ~ 30. While the best-£t values 


of th ese parameters are roughly consistent with the earlier analysis (e.g. iGofford et al 


20131 ). our model can further provide a geometrical and physical identihcation of the 


UFO in PG 1211-1-143 data rather than phenomenological interpretation as illustrated in 
Figure EJj. In fact, a recent discovery of the unambiguous P-Cygni-like prohle of the UFOs 
made by simultaneous X-ray observations with XMM-Newton and NuSTAR of a similar 
luminous quasar, PDS 456, indicates a simil ar spherically-exten ded wind geometry rather 
than a narrow collimated radial streamline (INardini et al.l 120151 ) in consistence with our 
MHD-driven view discussed for PG 1211-M43. 


Although in this paper we are focusing on the origin of the detected Fe Ka UFOs, our 
model winds extend over a large range in r, ^ and v. As such, they imply the presence of 
other charge states that contribute to the Fe Ka transition by including Fe xvill-Fe xxiv. 
We found that should these additional states be included in our analysis, the Fe Ka feature 
would have been much broader than seen in the data. Given our hts of Fig. 0^ and Table |2l 
one must surmise that the effective contribution to ls-2p transition from Fexvill through 
Fexxiv in the data ought to be very small (if not none). There are a number of remedies: 
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(i) It is conceivable that the intrinsically broad absorption featnre dne to ionized iron at 
all charge states conld be externally hlled by scattered resonant line photons which wonld 
snppress its otherwise broader signatnre. Any continnous wind model will inevitably come 
across this issne of the contribntion of states other than highly-ionized (e.g. H/He-like) 
ones, (ii) The radial wind density profile might be steep enongh to snppress the ionic 
colnmn at large distances. On the other hand, this solntion may not be consistent with 
the observed slow absorbers (i.e. warm absorbers) since they originate from large distances 
in this model, (hi) It is also probably that the fast absorbers (i.e. UFOs) conld be a 
collection of discrete (small) gas clonds along the LoS instead of a large-scale continuons 


flow fe.g. iMisawa et al.ll2014l) that mi ght also be in a constant pressn re equilibrium causing 


the suggested thermal instability (e.g. iHolczer. Behar fc Kaspil 120071) . While we note this 


long-standing question, this is beyond the scope of our current study. 

To compute the spectra of a truncated wind, we simply removed the contribution to 
the line feature by the self-similar section of the wind that originates aX R < ftRs- We also 
repeated the photoionization of a wind in which this section has already been removed, 
assuming always that the ionizing source is located at r = 0. This second calculation 
produced similar results with slightly smaller values for Mh, because of the slightly larger 
flux of ionizing radiation at the values of 6 considered. This could be ameliorated by a 
slight increase in the value of Uo- The constrained truncation radius Rt ~ 30i?5 in Table |2] 
is statistically favored in the context of our MHD-wind mod el particularly so as to su ppress 
the blue tail of the absorption feature. On the other hand, lOin.stini fc Proeal (1201211. for 


examp le, have considered a thermally-driven wind based on the model of iLuketic et ah 
fl2010l) . They found a relatively sharp blue edge of the line prohle without truncation due 


to non-monotonic prohles for wind streamlines and opacity along a LoS. This implies that 
a complex geometry of the wind also needs to be further explored by extending the model 
beyond self-similar limit. 

It is suggested from a long Suzaku observation that a similar fast X-ray absorber 
(by iron K-shell transition at an implied outfl ow velocity of n ^ 0.25c) in PDS 456 
exhibits rapid variability as short as ~ 1 week flGofford et al.l 1201411 . Their estimate of the 
absorber’s location in PDS 456 {r/Rs ~ 100 — 1800) is very similar to our estimate Tc 
in PG 1211-1-143. The current steady-state model is not appropriate for treating such a 
time-variability in its absorption features, but it is conceivable that the observed variable 
nature may be associated with the change in wind density (perhaps resulted from the 
change in mass-loading) and/or changing streamline configurations due to the variable 
magnetic helds. The model is in a good agreement with the data while there appears to 
be additional weak (intrinsic) absorbers at higher energies (~ 8 — 9 keV). These weak 
absorption structures can be due to the res onance series converging to the Fexxv edge (e.g. 
Kallman et al. 20041 : Tombesi et ah 2011al ). There can be also some contamination due to 
the presence of the background (instrumental) emission lines such as Gu Ka at 8 keV in the 
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EPIC-pn spectrum. 


In this paper we employ a well-studied semi-analytic wind model as a primary 
component. We feel that despite their simplicity, such models should not be dismissed 
off hand compared to large scale purely numerical simulations for a number of reasons: 
First, even the state-of-the-art simulations today have not yet provided the practical and 
direct observables addressed in this paper at an observationally-relevant level. Second, 
it is still computationally extremely challenging to include self-consistently multi-scale, 
multi-dimensional radiative transfer for plasma/atomic physics necessary to simulate the 
kind of transitions seen in UV/X-ray data while simultaneously covering a large spatial 
scale (say, ranging from 10 Schwarzschild radii all the way out to pc-scale) without 
suffering from numerical instability and boundary conditions susceptibility. In the next 
step a more self-consistent disk-wind morpholo gy needs to be considered by constructing 
a sop histicated (perhaps dynamical) model (e.g. Ohsuga et ah 20091: Ohsuga fc Mineshige 


2011 ) also incorporating detailed radiative transfer for spectral lines (e.g. 
2004 : Garcia et ah 2013 ). 


Kail man et ah 


In this preliminary calculation we assumed ma,o = 1 corresponding to the density of 
matter on the disk surface Uq = 5 x 10^^ cm“^. We have used slightly different values for Uo 
and noted their weak influence of the end results. We note that our assumed value is slightly 


higher than a hducial AG N value of 10 


10 3 


.e.g. 


Grenshaw et ah 2003: Tombesi et ah 


2011al: iGofford et ahl l2013l. and references therein) but its poss i ble range can be conside red 
as broad as 10^° < rio < 10^^ cm“^ (e.g. 


Laor & Netzer 

1989: 

George & Fabian 

1991 


Garcia et al.l 120131) and an accurate assessment requires a more realistic modeling of 
accretion disk physics and its response to photoionization process. 


The present analysis is based on a selected hducial wind structure that we have 
examined in our earlier work (F14). Within the three-parameter model spanned by 0, fcTbbb 
and Rt in this paper, we do not notice degeneracy in the best-ht spectrum. Gonsidering 
a complexity of magnetized disk-wind physics, however, it is conceivable that we may 
End another best-ht solution from slightly diherent wind conditions. Removing such 
potential degeneracy is in principle challenging since there is little a priori knowledge (at 
least observationally) of the underlying wind structure. Nonetheless, it will be possible 
to rule out some of the degenerate wind solutions by further including multiple ions of 
diherent charge states both at soft X-ray transitions (below 3-4 keV) and Fe-Ka transitions 
simultaneously since all these absorption signatures should be coupled in the context of our 
continuous disk-wind scenario. We thus plan to extend the current preliminary model to 
include the soft X-ray WAs to examine a coherent predictability of the model using those 
AGNs exhibiting both WAs and UFOs. 


We anticipate the upcoming missions such as Astro-H and Athena to contribute 
signihcantly to this goal by providing more detail on the Fe-K component of the wind as 
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well as soft X-ray absorbers, and thus to further clarify our picture of AGN structure. 

We are grateful to the anonymous referee for the constructive criticism to improve the 
quality of the manuscript. 
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Table 3. Various Characteristic Radii in This Model 


Symbol 

Description 

Ro 

The fiducial radius of the wind models 

Rs 

Schwarzschild radius 

Rt 

Disk inner truncation radius 

Tc 

Radius where a given ion becomes most opaque 

Risco 

Innermost Stable Circular Orbit 


All radii are defined in the cylindrical coordinates except for Vc which is measured along a LoS 
direction. 
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